Deducing spectroscopic factors from wave-function asymptotics 
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In a coupled-channel model, we explore the effects of the coupling between configurations on the 
radial behavior of the wave function and in particular on the spectroscopic factor (SF) and the 
asymptotic normalisation coefficient (ANC). We evaluate the extraction of a SF from the ratio of 
the ANC of the coupled-channel model and that of a single-particle approximation of the wave 
function. We perform this study within a core+n collective model, which includes two states of 
the core that connect by a rotational coupling. To get additional insights, we also use a simplified 
model that takes a delta function for the coupling potential. Calculations are performed for 11 Be. 
Fair agreement is obtained between the SF inferred from the single-particle approximation and the 
one obtained within the coupled-channel models. Significant discrepancies are observed only for 
large coupling strength and/or large admixture, i.e. small SF. This suggests that reliable SFs can be 
deduced from the wave-function asymptotics when the structure is dominated by one configuration, 
i.e. for large SF. 
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I. INTRODUCTION 

Historically, direct reactions have been an important 
source of information on nuclear structure. The develop- 
ment of radioactive-ion beams has rekindled the interest 
in those reactions by providing a unique way to study 
nuclei far from stability [![. Through direct reactions, it 
is indeed possible, with an appropriate reaction model, 
to extract quantitative information on the structure of 
nuclei bound for only milliseconds @]. For example, the 
orbital configuration of the valence nucleons in a num- 
ber of neutron-rich nuclei has been determined through 
transfer knockout [|| and breakup [f| reactions. 

The abrupt nature of direct reactions leads to lim- 
ited structural changes in the involved nuclei. Therefore, 
the cross sections for those reactions can be expressed in 
terms of overlap functions, which, strictly, involve overlap 
integrals of the full wave functions for the initial and final 
states of the nuclei @ . Information about overlap func- 
tions can be inferred from direct reactions. In particular, 
the square of its norm, called spectroscopic factor (SF), 
is supposed to be obtained from a comparison between 
theory and experiment. Some controversy exists about 
the use of SFs when discussing nuclear structure 0, [3 ■ 
Nevertheless, SFs are attractive in that they give an in- 
sight into the complexity of the nuclear states involved in 
the reaction. Moreover, they can be compared to results 
of electron inelastic scattering and of shell-model calcu- 
lations. Alternatively, it has been suggested that, since 



'Electronic address: picrrc.capcl@ccntralicns.net 
^Electronic address: daniclcwicz@nscl.msu.edu 
t Electronic address: nuncs@nscl.msu.edu 



direct reactions are mostly peripheral, they probe only 
the tail of the overlap function [jj, Gil • Accordingly, only 
the asymptotic normalization coefficient (ANC) of the 
overlap function should be extracted from data analysis. 

Currently there is no fully microscopic basis for calcu- 
lating overlap functions for nuclei with mass A > 10, let 
alone a reaction model that includes these microscopic 
overlap functions. Reaction data are usually analyzed 
within models where nucleons are grouped into few clus- 
ters, and in which overlap functions are substituted by 
wave functions obtained by solving a Schrodinger equa- 
tion in which the interactions between the clusters are 
simulated by local potentials. As the primary application 
of this work is the study of nuclear structure far from sta- 
bility, we consider the case of two-cluster systems, such 
as one-neutron halo nuclei, in which a neutron is loosely 
bound to the core of the nucleus. In that case, the overlap 
function is approximated by a single-particle wave func- 
tion, which describes the valence neutron evolving in a 
mean-field potential representing its interaction with the 
core. The cross section obtained with this single-particle 
approximation of the overlap function is rcnormalized to 
match experimental data and, from that rcnormalization, 
a SF or an ANC is deduced. In this work, we examine the 
validity of these assumptions by addressing the following 
questions: First, how accurately can one substitute the 
exact overlap function by a single-particle wave function? 
Second, if only the ANC can be extracted from measure- 
ments, can it be reliably related to the SF? 

Within a fully microscopic description of the nuclei, 
overlap functions satisfy a set of strongly coupled equa- 
tions. Formally, this set can always be reduced to one 
equation for the single function of interest, evolving under 
the influence of an effective single-particle Hamiltonian. 
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However, due to strong couplings between the different 
channels, that Hamiltonian is going to be highly non- 
local, in terms of its dependence on both position and 
energy. For instance, it will describe absorption of the 
probability flux from one channel into another channel. 
This transfered probability flux may be absorbed from 
one side of the nucleus and returned to the other side of 
the nucleus. Obviously, this involves a significant spatial 
nonlocality, up to the size of the nucleus. The associated 
energy dependence of the effective Hamiltonian may be 
tied to the time delay in returning the flux to the original 
channel. Needless to say that such type of effects can- 
not be properly accounted for in terms of a simple local 
potential. By contrast, coupling effects with exchange 
of probability flux between channels may be adequately 
simulated at a semiquantitative level in collective models. 

In this work, we consider the rotational-coupling model 
developed by Nunes et al. [THHH. ^ n m °del, the core 
of the nucleus is described as a deformed rotor, which can 
be in various excited states. This leads to a set of cou- 
pled equations for the core-nucleon wave functions sim- 
ilar to those obtained within microscopic models. The 
rotational model is applied to 11 Be, the archetypal one- 
neutron halo nucleus. This nucleus is described as a neu- 
tron loosely bound to a 10 Be core which can be in both 
its + ground state and 2 + excited state [llj. Note that 
similar models have been developed by other groups, see 
e.g. Refs. P. EH. 

In order to assess the validity of the single-particle 
approximation, and to understand the influence of the 
couplings upon an overlap function, its SF, and ANC, 
we carry out cross-comparisons within theory. We ana- 
lyze the more realistic rotational model in terms of the 
single-particle approximation. In particular we compute 
the ANCs within the rotational model and the single- 
particle approximation and we deduce a SF from their 
ratio. Confronting the deduced SF to that directly ob- 
tained within the rotational model, we can infer the reli- 
ability of the extraction of SFs from ANCs. This method 
could prove very valuable if ANCs can be efficiently mea- 
sured from direct reactions, like transfer 0] or breakup 
fioj . Note that besides the accuracy of the method, the 
extraction of SFs from ANCs is subject to other uncer- 
tainties. In particular, the geometry of the core-neutron 
potential is not known a priori. In this work, we focus on 
the validity of the method and disregard this latter uncer- 
tainty {9l. [l5L [iH] . Throughout our analysis, we always use 
the same potential geometry within both the rotational 
model and the corresponding single-particle approxima- 
tion. 

Besides the aforementioned rotational model [IJ [l2| , 
we also employ a simplified schematic model with two 
overlap functions generated with square-well potentials 
and coupled by a delta function. This simplified model 
may be largely solved analytically, which enables us to 
interpret qualitatively the results of the rotational model. 

An important issue not treated in the present work 
is the effect of short-range correlations, which can con- 



tribute to the lowering of SF. Those correlations are of- 
ten the focus of investigations in many-body theory of 
nuclear matter [13]. They generally affect the overlap 
function within the nuclear volume, just as do the cou- 
plings to collective states. Hence, although not exactly 
included here, short-range correlations will be qualita- 
tively simulated within the rotational model. 

Direct reactions are a useful tool to study nuclear 
structure far from stability. However, many uncertain- 
ties remain in the analysis of measurements. While sig- 
nificant advances have been made recently on the side 
of reaction theory [H|-[22j], the models employed in the 
analysis remain fairly schematic relative to the full many- 
body problem. On the structure side, the attention of 
microscopic models has been primarily directed at the 
treatment of short-range correlations. Traditionally, lit- 
tle attention was paid to the asymptotic part of the wave 
functions 0, [24|, which is essential in the analysis of 
direct reactions. More recently, there have been works 
focusing on the asymptotic behavior [25M27I] . However, 
the task of establishing this behavior may become very 
complicated for loosely-bound nuclei, and even more if 
multiple-channel configurations are involved. For these 
and other reasons, the single-particle model is likely to 
remain a work horse in data analysis in the near future. 
Therefore, assessing uncertainties associated with that 
model is important. 

In the next section we provide the necessary theoret- 
ical background for this study. We first introduce the 
notion of overlap functions and their couplings. We then 
present the various models used in this work: the single- 
particle model, the rotational model, and the schematic 
delta-coupling model. In Sec. IIII1 we provide quanti- 
tative results. First, we study the realistic conditions 
for the case of 11 Be. Second, we explore the limits of 
the parameter space in the effort to draw general con- 
clusions. Finally, we investigate these results within the 
more schematic delta-coupling model. Our work is sum- 
marized in Sec. IIVI 



II. THEORETICAL CONSIDERATIONS 
A. Overlap function, SF, and ANC 

Within direct-reaction theory, the cross sections for 
transfer, break-up, or knockout are usually expressed in 
terms of overlap functions @ . The reaction is indeed ex- 
pected to reveal the dominant configuration within the 
projectile structure, such as A — > B + n. Here and in the 
following, we consider a nucleus A that exhibits a strong 
two-cluster structure: a valence neutron n bound to a 
core B. To simplify matters, the symbols A and B will 
also stand for the mass numbers, i.e. B = A — 1. 

As an introductory step, we derive the equations satis- 
fied by the overlap wave function within a general many- 
body formalism. We will specifically refrain from dis- 
cussing the issues of spin, some antisymmetrization ef- 
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fects and center-of-mass corrections. With the exception 
of antisymmetrization, which cannot be considered ex- 
actly in a collective model, these will be accounted for in 
the practical calculations. 

Formally, the microscopic Hamiltonian describing the 
motion of the nucleons of nucleus A reads: 

A A 
i=l j>i=l 

where ti is the kinetic-energy operator for nucleon i and 
Vij describes the interaction between nucleons i and j. 
For simplicity, we do not mention three-body interac- 
tions. The states v A of A are the eigenstates of Ha, 

E$ A $ VA (r u ...,r A ) = H A * Va (n, . . . , r A ) , (2) 

where {t^} are the coordinates of the nucleons of A. Sim- 
ilarly, an (A — l)-body Hamiltonian, with eigenstates $„ B 
of energy , can be defined for the core B. As men- 
tioned earlier, the states of A may present a strong B-n 
cluster structure. In that case, it is worthwhile to de- 
scribe the relative motion between the core in its state 
vb and the valence neutron in terms of the overlap func- 
tion ip VA VB , which is nothing but the projection of fy VA 
onto the wave function describing the core state $„ B , 

^u AVB (r) = f dr! ■ ■ ■ dr A -i<f>* B (ri, . . . ,r A -i) 

x^ UA ( ri ,...,r A ^i,r), (3) 

where r is the coordinate of the valence neutron relative 
to the core. 

The spectroscopic factor is defined as the square of the 
norm of the overlap function: 

S VAVB = J dr\i> VAVB {r)\ 2 . (4) 

This SF represents the probability that, within the state 
v A of A, the neutron may be found in a combination with 
the core B in state vb- The SFs add up to unity over all 
the states of B, including those in the continuum, 

vb 

since that sum represents the square of the norm of ^! VA . 

To obtain the formal equations satisfied by the over- 
lap function, we project both sides of Eq. onto <& UB . 
Taking account of Eq. , we obtain 

K^VAVB = (^U B \H A \^ VA ) 

^(EvB+tr) i'u A ^B+Xi^fiB^XMB! ( 6 ) 
I 1 B 

where tr is the B-n kinetic-energy operator and 

B 

VvbUB = (^BlX^I^) (7) 

1=1 



are potential terms that couple different possible config- 
urations va vb within ^l VA . Let us now define the energy 
of the state v A relative to the B-n threshold (i.e. the 
negative of the separation energy) by 

e%=K- E o, (8) 
and the core excitation energy by 

<*b = E? B - E* • (9) 

In Eqs. ([SJ and ([5]), Eq is the ground-state energy of B, 
hence e®^ < for particle stable A and e„ B > 0. The set 
of equations ((HJ) can now be rewritten into the form 

{ £ v A ~ e v B ) VVa v b = (fa" "f" Kb v b ) VV>i v B 

+ V »b mb VVt mb ■ (10) 

Mb -/-v b 

The picture behind this set of equations is that of a neu- 
tron moving around B, which can be in its various states. 
The possible quantum numbers for the B-n relative mo- 
tion in each configuration v A vb are determined by con- 
servation laws and by the quantum numbers of the states 
va and vb- In particular, the orbital angular momentum 
/ of the B-n relative motion is determined by parity con- 
servation and angular-momentum couplings. 

The diagonal potential element V VB VB in Eq. (|10[) rep- 
resents the single-particle potential for n being in the field 
of B in its state vb- The motion of n is modified by the 
presence of other channels v A /j,b, which couple to v A vb 
through the non-diagonal potential elements V VB MB . In 
passing, we may note that, since the internucleon inter- 
actions vij are at most weakly nonlocal, both the single- 
particle and the coupling potentials (J7J) will be at most 
weakly nonlocal in r. In addition, this weak nonlocal- 
ity, combined with the short range of nuclear interaction, 
will lead to a form factor for the nuclear contribution 
to V VB VB that approximates the shape of the density of 
B. Outside the volume of B, the nuclear contributions 
to the potentials get suppressed and the wave equations 
simplify. For a B-n system, i.e. in absence of Coulomb, 
the radial part of the overlap function of a configuration 
with orbital angular momentum I exhibits the asymptotic 
behavior 

i>u A v B (r) > C VA VB i k ua u B h\ ' (i k va UB r) , (11) 

r— >oo 

where %K VAVB = ^2n Bn \e^ - e„ B |, with n Bn the B-n 

reduced mass, and hj is a spherical Bessel function of 
the third kind [28]. The normalization constant C VAVB 
appearing in Eq. (|11|) is the asymptotic normalization co- 
efficient. The function hj accounts for the effects of the 
centrifugal barrier and behaves as exp(— k Va Vb r) / n UA UB r 
for k va VB r 3> I- Hence, of all the different overlap func- 
tions, the one corresponding to the ground state of B 
(i.e. with e„ B = 0) dominates the asymptotic behavior of 
i\) VA VB , whereas the overlap functions corresponding to 
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high-lying states of the core barely stick out of the vol- 
ume of B. In the presence of couplings, any channel that 
satisfies the conservation laws could, in principle, con- 
tribute to the spectroscopic factor. However, in the case 
of loosely-bound systems, much of the wave function is 
expected to be outside the volume of B. The dominance 
of the B ground-state channel in the asymptotic region 
suggests that this channel also dominates the contribu- 
tions to the SF. 

Due to the couplings, the shape of the overlap func- 
tion tj) VA Ub differs from the solution of a single-particle 
Schrodinger equation in which the nuclear B-n interac- 
tion is simulated by a local potential tied to the density 
of B. This influence of the couplings upon the shape of 
the overlap function, and in particular on the connec- 
tion between SF and ANC, is the focus of the present 
work. As mentioned in the introduction, we do not solve 
the full many-body problem ([2]) but rather simulate the 
set of coupled equations ((TU)) using a collective model in 
which the core is described as a deformed rotor [ll[ . The 
results obtained within that rotational model will then 
be compared to those of a mere single-particle model in 
order to evaluate the sensitivity of the overlap function 
to the couplings. For a qualitative understanding of our 
results, we also study a more schematic model in which 
two channels are coupled through a delta interaction. 

For completeness, we lay out in Appendix [X] the for- 
mal reduction of the set (|10p to a single-particle equation 
with a nonlocal effective Hamiltonian. The nonlocality 
is dominated by couplings to low-lying excitations of the 
core, such as those investigated within the models here. 
That reduction can be, in fact, applied directly to the 
models we employ, and it is bound to yield the same 
shape of overlap functions we obtain. 

The different models considered here are presented in 
the following subsections. Each of them is particularized 
to 11 Be, the test case of our study. 



B. Single-particle approximation 

Within the single-particle (sp) approximation, only one 
configuration is considered. This comes down to neglect- 
ing the last terms in the r.h.s. of Eq. (|TU|) . which there- 
fore reduces to a mere one-body Schrodinger equation. 
After the factorization of the overlap function into its ra- 
dial and spin-angular y>ij parts, the single-particle 
Schrodinger equation reads 



the radial part of the B-n kinetic-energy operator reads 



(e 



Bn 
• 'A 



X j (r) = [i + V VBlj (r)}^ lj (r), (12) 



where n is the principal quantum number and j is the 
quantum number associated with the angular momen- 
tum obtained from the coupling of the orbital angular 
momentum I and the spin of the neutron. In Eq. (1121) . 



4 = {yij\t r \yij) 

h 2 / d 2 2d 

2/i Bn \ dr 2 r dr 



1(1 + 1) 



(13) 



As in the many-body case, the asymptotic behavior of 
the single-particle wave function, normalized to unity, is 



(14) 



where b n ij is the single-particle ANC. From Eqs. (fTTj) and 
(fT4"| . it is clear that the many-body overlap function is 
directly proportional to the single-particle wave function 
for large r. If a reaction probes only the ANC of the 
overlap function, and the assumption is made that the 
couplings to other states have little impact on the shape 
of the overlap function, then the SF deduced on the basis 
of the single-particle model amounts to 



QSp 



c v 



(15) 



This single-particle approximation is tested here in the 
case of 11 Bc. This typical one-neutron halo nucleus has 
two bound states: a l/2 + ground state and a l/2~ ex- 
cited state. If one assumes the 10 Be core to be in its + 
ground state, the corresponding single-particle orbitals 
are 2sl/2 and lpl/2, respectively. The nucleus n Be also 
exhibits a 5/2 + resonance, which is usually reproduced 
within the d5/2 partial wave coupled to the + ground 
state of the 10 Be core [2^|. If 10 Be can be in its 2 + excited 
state, then these orbitals couple to other configurations, 
as described in the following subsections. 



C. The rotational model 

We next consider a simple step beyond the single- 
particle approximation: we use a collective model where 
the core is deformed and allowed to excite ful4i~4|. The 
Hamiltonian of such a B-n system reads 



(16) 



where Hb is the internal Hamiltonian of the core, and 
now the effective interaction between the core and the 
neutron depends on the internal degrees of freedom of the 
core £. In the following, each level of the core considered 
in the calculation is identified by its spin / and parity 
tib- The corresponding wave functions and energies are 
denoted by &i* B and £/»b, respectively 



(17) 



The wave function of the B-n system is expanded in 
terms of these core eigenstates 



(18) 
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where the subscript i represents all possible quantum 
numbers that couple to the total angular momentum J 
and parity ir. For clarity, we omit the coupling coeffi- 
cients. In Eq. (fT5)l , the part of the wave function de- 
scribing the B-n relative motion, i.e. the equivalent of 
the overlap function is split into its radial and 
spin-angular 3fyj parts. Replacing Eqs. ([To]) and (fI5]l into 
the Schrodinger equation (J5|) , one arrives at the coupled- 
channel equations [TT| - |T3| : 



prise the following components 



0# 



e % W°\r) = [t l r + Vu(r)}rt°\r) 



(19) 



where the potential matrix elements Vij, responsible for 
the couplings between various components of the wave 
function, are defined by 

Va(r) = (m)yi(r)\VBn(r,0\yi(r)^(0)- (20) 

The coupled equations (fT9")l simulates those satisfied 
by the microscopic overlap function (JTOJ) . In the rotor 
model, the core is described as a deformed rotor fllj ]. 
Within this model, the internal coordinates of the core 
£ are the three Euler angles, and the functions &i* B are 
the Wigner rotation matrices D I MK , with M and K the 
projections of I in the lab and core intrinsic frames, re- 
spectively. For the application we have in mind, namely 
n Be, we consider only the first two states of the 10 Bc 
core: the + ground state, and the first 2 + excited state 
(e 2 + = 3.368 MeV). They are seen as the first two states 
of a rotational band of intrinsic projection K — 0. 

Consistent with a deformed core, the interaction be- 
tween the core and the neutron is described by a de- 
formed Woods-Saxon potential: 



Vws(r) = -Vi < 1 + cxp 



r - R(d, 



(21) 



in which the depth Vi may depend on the orbital angular 
momentum I. The radius R reads 



R(0,<f>) = Ro[l + [3Y 20 (e, </>)}, 



(22) 



where (3 characterizes the deformation of the core. Since 
this deformation is responsible for the couplings between 
configurations, j3 is also called coupling strength in the 
following. On top of this deformed central term, we also 
include the usual (undeformed) spin-orbit coupling term: 



Vso(r) = I ■ sVso 



ld_ 

r dr 



, r - R S o \ 
+ exp(— — j 



• (23) 



The radial components ip™^ are found by solving 
Eq. (|19j) . imposing bound-state boundary conditions and 
unit normalization for J* using the program face [T^j ■ 
Within this model, the wave functions of the l/2 + 
ground state and l/2~ first excited state of 11 Be Corn- 



er, rot 

W l/2+ 



= ci /2 (r)y il/2 (f)*o+(e) 

+ Tpl% /2 (r)y d5 /2(i-)®2 + (0 



*l/2" 



= <i/ 2 M3W 2 (r)$o+(£) 

+ <3/2M^3/ 2 (r)$ 2+ (0 

Asymptotically the radial behavior of is identical to 

nlj 



that of a single-particle wave function ip sp 



C 



nlj 



(24) 



but with a different ANC C^°* . Indeed, whereas ip^- are 
normalized to unity, tpffj have norms less than one due to 
their couplings with the other components. In addition, 
the radial dependence of ip n °ij ma y differ from its single- 
particle approximation because of these couplings. 
Within this model, SFs can be calculated directly 



rrot 
D nlj 



rotiVdr. 



0' 



(25) 



Comparing these "exact" SFs with approximation ([15")) 
provides a good test of the single-particle approximation 
used in most of direct-reaction models. 



D. The delta function coupling model 

To get a qualitative understanding of the results ob- 
tained in the rotational model, we have developed a more 
schematic collective model to describe the ground state 
of 11 Bc. Square wells are used as single-particle poten- 
tials and a delta function simulates the coupling term. 
For simplicity, the spin of the neutron is neglected. We 
consider two radial components il>o an d correspond- 
ing, respectively, to the valence neutron in the s wave 
coupled to the 10 Be core in its + ground state, and to 
the valence neutron in the d partial wave coupled to the 
first excited state 2 + of 10 Be. In this schematic model, 
the set of coupled equations (jTOj) reduces to 



e Bn y (r) = [t° r + V Bn (r)]4(r) 
+V s (r) xp 5 2 {r), 

( £ Bn -e 2 +)^(r) = [t 2 r + V Bn (r)]4(r) 

+Vs{r)$(r), (26) 

where Vbu is a square well of radius Rq and depth Vo, and 
Vg(r) = —(3 Vq Rq 5(r — Ro) with a coupling strength (3. 

Outside the potential, at r > Rq, the neutron wave 
functions exhibit the asymptotic form 



ip & {r) = Cq iK hy(iK r) 



(27) 
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where we have introduced the corresponding ANCs Cq 
and C 2 . Inside the potential well, at r < R , the wave 
functions are 

tf/ (r) = A j (k r) , 

Tp 5 2 (r) = A 2 j 2 (k 2 r), (28) 

where A a and A 2 are normalization constants, ji are 
spherical Bessel functions of the first kind (28j, HUq = 
yj2^ Bn {e Bn + V ), and hk 2 = yj2fjL Bn (e Bn - e 2 + +V Q ). 
At r = i?o, the overlap functions ipQ and i{j 2 are con- 
tinuous, but their derivatives are not, due to the delta 
coupling. This discontinuity is proportional to the value 
of the wave function in the other channel. 

The conditions at r — Rq may be combined into one 
equation relating e Bn , Vq and /?: 

{[\og{r^)]'\ Ro+ ~[\o g {r^)\\ Ro _} 

xjflogCr^'l^-ClogCr^)]'!^} 

= ^ Bn ^RoV y (2g) 

Here, Rq- and Ra+ stand for limits taken from below 
and above Rq. At the extreme of a vanishing coupling, 
P = 0, the equation (|2U|) will be satisfied when either of 
the l.h.s. factors vanishes. Vanishing of cither of those 
factors is indeed the standard square- well condition for a 
bound state, when no coupling exists. 

In the presence of a coupling, we find it easiest to fix 
Vo, which sets the value for the l.h.s. of (|29[) . and read 
off p from the r.h.s. Thereafter, we can determine the 
ratios of the constants Ai and C/ and, finally, absolute 
values for those constants from normalization of the wave 
function to unity. 

Besides the fact that it can be solved analytically, the 
delta-coupling model presents the advantage that it en- 
ables one to easily interpret the features of the set of cou- 
pled equations (|26|) in terms of flux of probability trans- 
fered from one channel to the other. As the wave for a 
given / encounters the boundary of the potential gener- 
ated by the core, angular momentum may be exchanged 
with the core and flux of probability can be transferred 
between the channels. Of course the net probability flux 
is conserved within the whole wave function, but not in- 
side a single channel: When starting from a solution for 
1 = 0, without coupling, and then enhancing the cou- 
pling, the picture is that of the coupling acting as an 
antenna, or source term, that radiates into the I = 2 
channel. For a wave traveling onto the coupling, the flux 
from the s channel to the d one is 

^K(^k_-<k)} 

= R 2 ^PR V $ WVo} k. ( 3 °) 

which shows that the discontinuity in the derivative of the 
wave function due to the delta coupling is related to the 



transferred flux of probability. In the stationary states 
we are interested in, the flux leaking out of a channel is 
exactly balanced by the flux coming in. 

The SFs are the contributions from the two compo- 
nents to the square of the norm of the wave function, 

/>oo 

Sf = / \^f\ 2 r 2 dr. (31) 
Jo 

Again, we can compare those SFs coming directly from 
the solution of Eq. f|26[) with their single-particle approx- 
imation (|15[) . The change in the shape of the wave func- 
tion, compared to the single-particle approximation, is 
due to the fact that the derivative is discontinuous at 
r = Rq. From the set ([26|) . it follows that the magnitude 
of the I = 2 wave function is linear in the coupling P, for 
small p. With this, the discontinuity in the derivative of 
I = wave function is quadratic in p. In consequence, 
any discrepancy between the exact spectroscopic factor 
and that estimated in the single-particle approximation 
should be quadratic in the coupling. At the same time, 
deviation of the spectroscopic factor from unity, for 1 = 0, 
should be quadratic in the coupling as well. In the end, 
the discrepancy in the spectroscopic factor is expected to 
be linear in the deviation of the dominant spectroscopic 
factor from unity. 



III. RESULTS 
A. Realistic n Be 

There are three well known states in n Be to which we 
fit the B + n interaction: the l/2 + ground state (g.s.) 
with an energy e?^+ = —504 keV, the l/2~ first ex- 
cited state (e.s.) with ef£_ = -184 keV, and the 5/2+ 
resonance located at = 1.274 MeV in the contin- 

uum. Our starting point for the undeformed interac- 
tion (i.e. the uncoupled case) is the potential developed 
in Ref. [29[. The radius and diffuseness of the Woods- 
Saxon form factor are chosen equal to Rq = 2.585 fm, and 
a = 0.6 fm, respectively. The spin-orbit depth is fixed at 
Vso = 21.0 MeVfm 2 , and its radius at R SO = 2.585 fm. 

We wish to explore the effect of the single-particle ap- 
proximation as a function of the coupling strength. Thus, 
although the quadrupole deformation P has a definite 
value (P = 0.67), which is related to the B(E2; 0+ -> 2+) 
of 10 Be [ll[, we use it as a parameter. For each P, the 
depths of the potential in the s wave (V s ) and d wave 
(Vd) are adjusted to reproduce the energy levels of the 
two positive-parity states. The depth is chosen the same 
in all the odd waves (V p ) and is adjusted to the energy of 
the e.s. The radius Rq is modified as a function of P to 
conserve the volume of the core. The spin-orbit coupling 
term is kept constant. 

In FigfT] we show the spectroscopic factors S^°* ob- 
tained from the solution of the coupled-channel equations 
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FIG. 1: (color online) Spectroscopic factors S^fc (|25|l and 
single-particle approximation $jjy (|15p as a function of de- 
formation. Note that the diamonds are superimposed on the 
circles. 



as a function of deformation. Results for the g.s. (cir- 
cles) and the e.s. (up triangles) are both included. The 
effect of deformation is very different in these two states 
although they are both loosely bound, and both corre- 
spond to a mixing of three configurations. The S^\, 2 
in the g.s. remains above 80% even at large deforma- 
tion. The remaining strength is shared between both d 
configurations. On the contrary the 5^°* , 2 in the e.s. 
drops rapidly with deformation down to 50-60%. The 
remaining part of the strength is exclusively in the p3/2 
configuration. No strength is found in the /5/2 config- 
uration. These results are in agreement with those of 
Refs. [U, HH, EH- For both states, there is no config- 
uration crossing: only one configuration dominates for 
all deformations. In both cases, the dependence of SFs 
on /3 is strongly non-linear. Interestingly, at small /?, and 
when there is little admixture, 1 — S^\, 2 behaves roughly 
quadratically in j3, as discussed at the end of Sec. Ill Dl 

We compare these SFs with the values that would have 
been obtained under the single-particle approximation 
(|15jl . For the g.s., this approximation (diamonds) works 
very well: no matter how large the coupling strength, 
^2si/2 ~ ^2°i/2- F° r the e.s. (down triangles), only the 
general trend of * s reproduced by S^p 1 / 2 - Although 

the single-particle approximation is valid at small defor- 
mation significant discrepancies appears for /? > 0.2. 

These results can be further illustrated by plotting the 
2s 1/2 and lpl/2 components of the wave function for 
various deformation parameters (Figs. [5] and [3]). These 
radial functions have been normalized to one to make the 
comparison easier. We see that the couplings simply re- 
arrange the interior of the sl/2 contribution: increasing 
j3 moves strength from the second peak to the first. How- 
ever the asymptotic part (beyond 5 fm) is left totally un- 
changed, explaining the validity of the single-particle ap- 
proximation (fT5j) . For this loosely-bound l/2 + state, we 
find that 50% of the spectroscopic strength comes from 
radii smaller than the interaction radius. This illustrates 
that even for loosely-bound systems, the internal part 



0.4 



0.3 



0.2 



0.1 



1 1 1 


1 1 1 1 


['' X 

//' \ 


P = 0.2 


✓V B \ 


P = 0.4 


&*. § \ 


= 0.6 " 


1 V P 


Nv p = 0.8 


/ F (? 

U i i 


i i i i 



6 8 
r (fm) 



10 12 



14 



FIG. 2: (color online) Radial part of 2sl/2 overlap functions 
normalized to unity. 
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FIG. 3: (color online) Radial part of lpl/2 overlap functions 
normalized to unity. 



of the overlap function significantly contributes to the 
SF. The validity of the single-particle approximation is 
therefore surprising. Although couplings significantly af- 
fect the interior of the overlap function, which contributes 
to a large part of the SF, the approximation (fT5"j) . based 
solely on asymptotic characteristics of the wave function, 
gives a precise estimate of the SF. 

As opposed to the 1/2+ state, the couplings seem to 
affect the wave function of the l/2~ state all the way 
out to large distances (see Fig. [3]). In this case, strength 
is moved outwards when j3 increases. This explains the 
deviation between Sl°\^ 2 and its single-particle approxi- 
mation S'lpj/a observed in Fig. [T] 

These results obtained for a realistic description of 
n Be are rather surprising. First, it was expected that 
at some point, e.g. for large deformation, the single- 
particle approximation would break down. However, as 
shown in Fig. [1] the single-particle approximation seems 
to be valid even at large /3, especially for the ground 
state. Second, large differences are observed between the 
two states. For the ground state, the 2sl/2 configuration 
dominates and the single-particle approximation is nearly 
exact. On the contrary, the excited state presents a large 
admixture of the lpl/2 and lp3/2 configurations, and 
the single-particle approximation of its SF is less precise. 

The difference in the validity of the single-particle ap- 
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proximation (|15|) can be understood in terms of the dif- 
ferent admixtures that arc obtained in the g.s. and the 
e.s. This suggests that the single-particle approximation 
is less valid at large admixture, i.e. when the correspond- 
ing SF is small. This agrees, at least qualitatively, with 
the arguments expounded in the description of the delta- 
coupling model. 

One possible explanation for the difference in the con- 
figuration admixture between both states is the presence 
of a bound state in the p3 /2 p artial wave generated by 
the core- neutron potential [29(. The existence of such a 
state could attract probability flux to that configuration 
and so enhance the admixture in the e.s. Accordingly, the 
low admixture observed in the g.s. would be explained by 
the absence of bound state in the d waves. To better un- 
derstand these results, we perform similar calculations 
extending the range of parameters to extreme values. 



B. Extreme couplings 

To study the difference in admixture observed between 
the two bound states, we perform calculations for a hy- 
pothetical g.s. of u Be with different values of Vd to see 
if and when the presence of a bound state in the d wave 
affects the admixture in that state. To do so, we repeat 
the analysis performed in Sec. IIII Al for the g.s., but now, 
instead of adjusting Vd to the resonant state, we use it 
as a free parameter and adjust V s to reproduce the cor- 
rect g.s. energy. For simplicity, the spin of the neutron 
is neglected. Results are shown in Fig. for deforma- 
tion parameters (3 = 0.1 (circles), 0.4 (squares), and 0.8 
(diamonds). The most striking feature is the pronounced 
drop of S^f at V d « 80 McV. At this depth, the single- 
particle potential hosts a d bound state with an energy 
that corresponds approximately to sfj^ + — e 2 + ■ In that 
condition, we have a nearly degenerate system, where the 
g.s. can be simultaneously in both 2s and Id configura- 
tions. It is only the presence of this deep d bound state 
that explains the large admixture. Indeed at lower or 
larger Vd, the admixture vanishes and > 0.8. This 
behavior holds for all coupling strengths, though the re- 
gion of large admixture increases with deformation j3. 
Note also the discontinuity of our results, especially at 
large /3. In the large-admixture range, parameters can- 
not be adjusted to reproduce the binding energy of the 
system due to numerical limitations. 

This result confirms the hypothesis formulated at the 
end of Sec. IIII Al It shows that large admixture is to be 
expected near degeneracy, i.e. when the coupled wave 
hosts a bound state close to the energy £j" — e/» B . 
However, if that bound state is located above or below 
that energy, no significant admixture should be observed. 
This explains the small admixture observed in the g.s. of 
the realistic 11 Be. In that case Vd is adjusted to repro- 
duce the 5/2 + resonance. In the single-particle potential, 
all d states are thus in the continuum, far from the en- 
ergy at which there would be degeneracy. The admixture 
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FIG. 4: (color online) (a) Spectroscopic factors as a function 
of Vd- (b) Relative effect of the single-particle approximation 
as a function of Vd- (a = 0.6 fm) 



observed in the e.s. of the realistic n Be is also explained 
by this effect. The single-particle parameters indeed lead 
fortuitously to a lp3/2 state bound by about 3 MeV [29j ] . 
which corresponds approximately to sfj^- — £2+ > i-e. the 
energy at which there is degeneracy. Note that the large 
admixture observed in the e.s. of 11 Be is unphysical as 
the p3/2 orbital is Pauli blocked to the valence neutron 
[30j | . Because of the large centrifugal barrier, the 1/5/2 
state is located high in the continuum, and so does not 
couple to the p configurations. 

Besides aiding to understand the results obtained in 
Sec. IIII Al these manipulations enable us to induce large 
admixture in our two-cluster rotational model. In this 
way we can study the validity of the single-particle 
approximation in nuclei with large fragmentation of 
strength that usually require more advanced structure 
models. In Fig. H£b), the ratio of the single-particle es- 
timate S S 2 P S (USD to its "exact" value S T 2 f |(2SJ), is plotted 
as a function of Vd- The approximation works perfectly 
at small coupling strength and its validity decreases as 
P increases. However, no matter how large the coupling 
strength, the single-particle approximation remains valid 
outside the large-admixture range with less than 10% er- 
ror. This confirms that the single-particle approximation 
breaks down mostly at large admixture. 

Fig. [5] shows the equivalent of Fig. |4] but for a much 
smaller diffuseness a = 0.1 fm. We perform this set 
of calculation not only to test the validity of our con- 
clusions in extreme cases, but also to ease the com- 
parison between this rotational model and the delta- 
coupling model, which corresponds to a nil diffuseness 
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(see Sec. IIIIC) . In this case, the large-admixture range is 
increased. Therefore the breakdown of the single-particle 
approximation happens sooner and to a larger extent. 
This is due to the smaller diffuseness, which produces a 
more abrupt change in the radial behavior of the wave 
function at the surface and so enhances the couplings 
between the various configurations. 

To summarize this analysis, we plot in Fig. [6] the 
single-particle approximation against the "exact" 
S 1 ^* for all the cases explored in this section. These 
correspond to the deformations (3 = 0.1, 0.4, and 0.8 
(circles, squares, and diamonds, respectively) and two 
diffuscnesses a — 0.1, and 0.6 fm (open and solid sym- 
bols). The dashed line corresponds to = and the 
deviation from this line estimates the error of the single- 
particle approximation. The main information conveyed 
by Fig. [6] is the general agreement between the "exact" 
SF and its single-particle approximation. Although there 
are some deviations from the = line, a large 
correctly predicts a large 5*2°*, while a small is 
usually obtained for small S^jf ■ Besides this qualitative 
information, Fig|B]emphasizes the fact that for small cou- 
pling strengths, the agreement between and is 
excellent, even for large admixture with other configura- 
tions, i.e. for small SF. It also shows very clearly that the 
deviation becomes larger as the coupling strength j3 in- 
creases. As noted earlier, this effect is more pronounced 
for smaller diffuseness. 

The largest discrepancy between the prediction S^ s 
and the "exact" 5*2°' is observed for small SF, where the 
single-particle approximation tends to significantly un- 
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derestimate the SF. To understand this effect, we display 
in Fig. [7] the radial component of some 2s overlap func- 
tions obtained in the most unfavorable case, i.e. /3 = 0.8 
and a = 0.1 fm. These functions are labeled by the value 
of Vd they correspond to (cf. Fig. [5]). The single-particle 
wave function (i.e. /? = 0) is shown as well for compari- 
son (thick full line). As already observed in Fig. [21 the 
couplings between the configurations affects the overlap 
function. In most of the cases these changes are lim- 
ited to the interior of the overlap function, leaving the 
asymptotics (nearly) unchanged. However, for some ex- 
treme cases, they extend much beyond the range of the 
single-particle potential. In particular, the node of the 
overlap function may be pushed so much outwards that 
the ANC Clf nearly vanishes (see, e.g., V d = 72 MeV). 
This leads to the very small estimates observed in 
Figs. [3|b), [3b), and[H whereas there remains some sig- 
nificant probability strength in the interior of the wave 
function, i.e., a non-zero ■ These rapid variations in 
the wave function also explain the difficulty in adjusting 
the parameters of the model within the large-admixture 
range (see Figs. 0] and [5]). However, these extreme cases 
seem to occur only for a combination of large coupling 
strengths and significant admixture. 



C. Analysis with the S coupling 

To expand our perspectives on the results obtained 
within the rotational model, we perform a similar analy- 
sis using the delta-coupling model developed in Sec. Ill Dl 
We consider here the g.s. of a 11 Be-like system. We take 
into account the main component 2s with the core in its 
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FIG. 7: (color online) Radial part of the s overlap functions 
normalized to unity obtained for /3 — 0.8 and a — 0.1 fm. 



FIG. 8: Spectroscopic factors as a function of deformation for 
the delta-coupling model. 



+ ground state coupled to the Id component with the 
core in its 2 + excited state. In a hrst step, the same 
square-well potential is taken in both s and d partial 
waves. 

In Fig. [5] we show a plot equivalent to Fig. []] for the 
delta-coupling model. We find that Sq (solid line) de- 
creases with the coupling strength j3 to very small val- 
ues. Note that the quadratic dependence of 1 — Sq in /3 
is only observed for very small coupling strength, namely 
(3 < 0.02. For f3 > 0.4 the dominant component in the 
1/2+ state is no longer the 2s component but the Id 
one. This larger admixture as compared to the rotational 
model may be due to the node of the 2s. In the rotational 
model, the combination of that node and the finite ex- 
tension of the coupling potential might sufficiently reduce 
the source term in the d-w&ve equation to limit the ad- 
mixture. On the contrary, the nil extension of the delta 
coupling avoids this cancellation effect, which could ex- 
plain the larger admixture. Another explanation is the 
intrinsically larger coupling strength in the delta cou- 
pling due to its nil diffuseness. We have indeed seen in 
Sec. lIII Bl that decreasing the diffuseness increases the ad- 
mixture between configurations (compare Figs.^a) and 
[SJa)). Whatever the reason for that larger admixture, 
it leads to a less reliable single-particle approximation 
(dashed line), as already seen in the rotational model. 
The Sq P now deviates significantly from the Sq (the er- 
ror can be as large as a factor of 2). Nevertheless, it 
reproduces qualitatively the general trend of the "exact" 
SF. Moreover, at small admixtures (i.e. for Sq > 80%), 
the single-particle approximation is in perfect agreement 
with Sq. 

To confirm the degeneracy effect observed in the ro- 
tational model, we repeat the analysis performed in 
Sec. IIIIBI The depth of the square well in the d wave 
Vd is varied while the depth in the s wave is adjusted to 
fit the 10 Be-n binding energy. The results of this analysis 
are displayed in Fig. |H1 where Sq is plotted as a function 
of Vd for various coupling strengths f3 = 0.05 (full line), 
0.2 (dashed line), and 0.4 (dotted line). These results are 
very similar to those obtained for the rotational model, 
and in particular for the small diffuseness (see Fig.[3Ja)). 
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FIG. 9: (color online) Spectroscopic factors as a function of 
Vd for the delta-coupling model. 

Little admixture is observed outside a range centered on 
Vd ~ 75 MeV, at which value the single-particle potential 
hosts a d bound state at e Bn — e 2 + ■ As in the rotational 
model, the width of that large-admixture range increases 
with the coupling strength. We also observe disconti- 
nuities in the SF similar to those observed in Fig. [5ja). 
These results confirm the degeneracy effect observed in 
Sec. IIII Bt large admixture takes place when the poten- 
tial well in the coupled channel is quantum-mechanically 
fit to host a bound state at the right energy. This effect 
seems therefore general and not due to some artefact of 
the rotational model. 

In Fig.rjDl the single-particle approximation Sq P (fT5)) is 
compared to the "exact" Sq (fBTj) . This figure is very sim- 
ilar to Fig. [6] It confirms that the single-particle approx- 
imation predicts a SF in qualitative agreement with the 
"exact" one. The prediction is more reliable for small ad- 
mixture, i.e. large SF, as already seen in Fig. [6] At these 
values, the discrepancy between the "exact" SF and its 
single-particle approximation roughly follows 

S s - S S P cx 1 - S 5 , (32) 

in agreement with the reasoning of Sec. Ill Dt This er- 
ror increases with the coupling strength ft. Note that a 
similar, though less obvious, result is observed in Fig. [51 
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FIG. 10: (color online) Comparison of the single-particle ap- 
proximation Sq P (|15[) with the SF obtained within the delta- 
coupling model So (|31[) . 



IV. DISCUSSION AND CONCLUSIONS 

Spectroscopic information about exotic nuclear struc- 
tures, like halo nuclei, can be inferred from direct reac- 
tions 0-13. Usually a SF and/or an ANC are extracted 
from the analysis of these experiments. Most of these 
analyses are performed within a single-particle approxi- 
mation of the projectile wave function. In order to eval- 
uate the validity of this approximation, we have studied 
the influence of coupling between configurations upon the 
overlap wave function. Keeping in mind the extraction of 
SF from peripheral reactions (3, [l(| , we have focused our 
study on the effects of these couplings on the SF and the 
ANC. Assuming that an ANC can be reliably extracted 
from direct reactions, we have analyzed how accurately 
a SF can be deduced from it. 

In this work, we have considered a collective model 
of the nucleus in which a valence neutron is bound to a 
deformed core allowed to excite. In particular, we have 
used the rotational-coupling model of Nuncs et al. [O, 
H3 |. in which the core is described as a deformed rotor. 
In this model the various states of the core are assumed to 
be part of a rotational band. A qualitative delta-coupling 
model has also been used to investigate the results of the 
rotational model. Calculations have been performed in 
the case of the typical one-neutron halo nucleus 11 Bc. 

Interestingly, this analysis shows that even a small 
coupling may affect significantly the overlap function. 
However, the changes compared to the single-particle 
wave function appear mostly in the interior, leaving the 
asymptotics nearly unchanged (see Fig. [5]). This sur- 
prising result suggests that probing only the tail of the 
wave function may give a reasonable estimate of the SF, 



contrarily to what has been assumed in Ref. [l(|. To 
explore this idea, we vary the deformation of the core 
that acts as a coupling strength. The single-particle ap- 
proximation (fT5)) is then compared to the "exact" SF 
obtained within the rotational model. This approxima- 
tion is obtained from the ratio between the ANC of the 
rotational-coupling model and that of the single-particle 
model. The former acts as the actual ANC of the system, 
supposedly measured from peripheral reactions. The lat- 
ter corresponds to the one obtained with a single-particle 
description of the nucleus, which is used in most of the 
reaction models. The prescription (fT5)) therefore sim- 
ulates the procedure performed in extracting SFs from 
experimental data. 

The results of our analysis show that the single-particle 
approximation often provides a reliable estimate of the 
SF. This estimate is at least in qualitative agreement 
with the coupled-channel models: It predicts large (re- 
spectively small) SFs, when large (respectively small) SFs 
are obtained from the coupled-channel models. For very 
small coupling strength, this estimate is very accurate. 
On the contrary, for large coupling strength, and in par- 
ticular for large admixture of the configurations, i.e. small 
SF, the single-particle approximation is much less reli- 
able. 

Being obtained within two independent structure mod- 
els, this result seems quite general. To understand it on 
firmer grounds, let us look back at the coupled equa- 
tions satisfied by the overlap function (|10[) and compare 
them with the single-particle equation (|12[) . The overlap 
function i\> VAVB will depart from its single-particle ap- 
proximation only when the coupling terms in Eq. (|10[) 
become significant, i.e. when at least one of the products 
V VBflB ip VAflB is large. This will happen for large coupling 
strengths, and/or when the population of the vaI^b con- 
figuration is large in the range of V UBflB , i.e. when there 
is a significant admixture between the various channels. 
In these two cases, the coupling affects the overlap func- 
tion so much that the single-particle approximation no 
longer holds. This reasoning also illustrates the possible 
influence of a node in one of the overlap function upon 
the admixture. Such a node may indeed cancel out the 
effect of the coupling if it is located within the range of 
Vvbub- I R the future, we plan to further investigate this 
effect. 

This study suggests the simple rule of thumb that large 
SFs deduced from ANCs are reliable. On the contrary, 
single-particle analysis of data estimating small SF, i.e. a 
large admixture between configurations, should not be re- 
lied on. Note that in any case some uncertainty remains, 
and that even large SF estimates can be off by 10-20%. 
Therefore, these estimates must be taken for what they 
are — estimates — and not precise, unquestionable values. 

The accuracy of the single-particle approximation de- 
pends also on the geometry of the potential chosen to sim- 
ulate the mean field of the core. In practice this geometry 
is unknown. We have not considered this uncertainty in 
the present analysis. Fortunately, as observed by Sparcn- 
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berg et al. [3l[ , the single- particle ANC for loosely-bound 
states is not strongly dependent on the potential geom- 
etry. In particular for an s valence neutron, it can be 
efficiently estimated from only the binding energy of the 
system (see Eq. (20) of Ref. (3l|). 

In our study, we consider a collective model to simu- 
late microscopic effects. To evaluate the sensitivity of our 
conclusions to this approximation, it would be interesting 
to repeat this analysis with microscopic models, such as 
the microscopic cluster model [H], the fermionic molec- 
ular dynamics [33l or the Green's function Monte 
Carlo model [2j], which all can correctly describe the 
asympotics of overlap wave functions for loosely-bound 
nuclear systems. 
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Appendix A: Effective Hamiltonian in a 
single-particle equation 

Given the linearity of the set of equations (|10[) for the 
overlap functions, any overlap function is linear in any 
other of those functions. Any selected equation may be 
then formally closed by expressing other overlap func- 
tions in terms of the function of interest. The cost is 
in a potentially strong nonlocality, both in position and 
energy, for the effective Hamiltonian within the single- 
particle equation for the chosen overlap function. Also, 
it may be difficult to construct the effective Hamiltonian 
in practice. 

Let the chosen overlap function correspond to the state 
vb of the core B. First, energy eigenvectors need to be 
found in the remaining space of overlap functions, in the 
absence of coupling to vb- With different energy values 
and eigenvectors denoted with index a, an eigenvector 
rf B a solves the set for energy E VB a , 

I rpU B a \ V B a _ f V B a . \ ^ ir V B a 

(for fi B + v B ) . (Al) 

Note that, without the coupling potentials V^ BTB , 
the states {rf Ba } would just represent energy eigen- 
states within the single-particle potential, coming in 
near-degenerate sets, because of the expected weak- 
dependence of the potential on fis ■ Given the coupling to 



i>v A v B of the overlap functions tp VA fi B i where fis ^ vb, 
the set of the latter functions can be next expressed as a 
combination of the vectors {rf Ba }. Upon inserting the 
expressions for the functions ip„ A ^ B into the equation for 
the overlap function il) VAUB , we arrive at the closed equa- 
tion: 

(e£T ~ VW B (r) = (tr + V„ B VB (r)) i> VAVB (r) 

+ J dr'AV VBVB (rS;e%)4>w B (r'), 

(A2) 

where the modification of the single-particle potential is 
AV UBVB (r,r';E)= £ V VBftB (r) 

V- B =£v B 
l B =£v B 

xGll 7B (r,r';E)V JBVB (r'), (A3) 

and the Green's propagator within the space complemen- 
tary to vb is 

a 

x<rW [<| a (0]* • (A4) 

The results are written here assuming local 2-body inter- 
actions. 

The nature of the contributions to the Green's func- 
tion, and to the modification of the potential AV, from 
different excitations of the core, will depend on the core 
excitation energy. High-lying excitations will yield con- 
tributions where propagation both to the outside and 
inside of the nucleus represents a tunneling. Those con- 
tributions will not yield much nonlocality for AV, nei- 
ther as a function of position nor energy. On the other 
hand, for low-lying excitations, the propagation within 
the nuclear volume may be largely uninhibited, while still 
turning to tunneling in the exterior. The latter excita- 
tions will contribute nonlocalities to AV, both in terms 
of energy-dependence and position, with latter nonlocal- 
ities possibly extending across the nucleus. The differ- 
ent limits of nonlocality may be, in particular, easily as- 
sessed in the delta-coupling model. In this paper, the 
focus is on the impact of coupling to low-lying states, on 
the overlap functions and on the deduced SFs. We may 
mention, though, that many microscopic interactions Vij 
with strong short-range repulsion may yield unphysical 
single-particle potential V VB VB in absence of any modi- 
fication in the form of AV. For such interactions, some 
level of renormalization of the potentials V VB ^ B , for low- 
lying states, through coupling to high-lying states may be 
necessary to start with. The consequence in the possible 
ill-definition of SFs, due to the latter type of couplings, 
has been indicated in Ref. @- 
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